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Abstract 

We consider a system in which a classical oscillator is interacting 
with a purely quantum mechanical oscillator, described by the La- 
grangian L = ^x 2 + ^A 2 — ^(m 2 + e 2 A 2 )x 2 , where A is a classical 
variable and x is a quantum operator. With (x(t)) = 0, the relevant 
variable for the quantum oscillator is (x(t)x(t)) = G(t). The classi- 
cal Hamiltonian dynamics governing the variables A(t), Ha (t), G{t) 
and n"G*(t) is chaotic so that the results of making measurements on 
the quantum system at later times are sensitive to initial conditions. 
This system arises as the zero momentum part of the problem of pair 
production of charged scalar particles by a strong external electric 
field. 



1 Introduction 

The definition and observation of chaotic behavior in classical systems is fa- 
miliar and well understood [I]]. However the proper definition of chaos for 
quantum systems and its experimental manifestations are still unclear ||]. 
Here we present a simple model of a coupled quantum-classical system and 
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introduce a new phenomenon that we call semi- quantum chaos. In a classical 
chaotic system such as the weather we are accustomed to situations where 
there is lack of long time forcasting because of the sensitivity of the system to 
initial conditions. The simple model we present here has the unusual feature 
that one has to give up long term forcasting even for the quantum mechanical 
probabilities, as exemplified by the average number of quanta at later times. 
The complete dynamics of the coupled quantum and classical oscillators is 
described by a classical effective Hamiltonian that is the expectation value 
of the quantum Hamiltonian. This effective Hamiltonian displays chaotic 
behavior, and thus the parameters that describe the quantum mechanical 
wave function (and hence expectation values) are sensitive to initial condi- 
tions. Chaos in dynamical systems with both quantum and classical degrees 
of freedom has been noted in more complicated systems and in a different 
context by other authors, (see e.g. M). 

We consider a system in which a classical oscillator is interacting with a 
purely quantum mechanical oscillator described by the Lagrangian, 

L = -x 2 + -A 2 - -(m 2 + e 2 A 2 )x 2 , (1) 

with equations of motion given by 

x + (m 2 + e 2 A 2 )x = (2) 

A + e 2 x 2 A = 0. (3) 

The Hamiltonian is 

H= 1 -p 2 + 1 -Il 2 A + 1 -(m 2 + e 2 A 2 )x 2 , (4) 

where pit) = x(t) and = A(t). We take x(t) to be a quantum operator and 
A(t) to be the amplitude of the classical oscillator. We require [x(t),p(t)} = i. 
We now introduce time-independent Heisenberg representation creation and 
destruction operators, a and a* , by the Ansatz 

x(t)=f(t)a + t(t)a^ (5) 

and we note that if f(t) satisfies the Wronskian condition 

iT(t)f(t)-r(t)f(t)] = i, (6) 
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then a and a) satisfies the relation [a, a'] = 1. From (0) and (||), we find that 
fit) satisfies the equation of motion 

f + (m 2 + e 2 A 2 )f = 0, (7) 

with the normalization fixed by the Wronskian condition (||). We can satisfy 
these two equations by the substitution 



f(t) = exp 



tt{t')dt' 



where Q(t) satisfies the nonlinear differential equation 
with 

u 2 (t) = m 2 + e 2 A 2 (t) . (9) 

Now, we choose the initial state vector at t = to be the ground state 
of the operator n = a^a, |^(0)) = |0), where a\0) = 0. Then, from (^), the 
average (classical) value of x(t) and pit) is for all time, (x(t)) = and 
(p(t)) = 0. However, the quantum fluctuations of x(t) are non-zero and are 
given by the variable G(t), 

G(t) = (At)} = \f(t)\ 2 = ^ f y (10) 



Then, from (^), it is easy to show that G(t) satisfies 

lu 2 = 0. (11) 



1 fG\ 1 fd\ 2 1 



2\GJ A\GJ AG 2 
In addition, we find that 

The expectation value of Eq. (ffl) becomes a new effective Hamiltonian 

#eff = (H^)) 

= ^ + 2U 2 G G + ^ + 1 -(m 2 + e 2 A 2 )G. (13) 
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The conjugate momenta are 

n G = ^, n A = A, (14) 

This classical Hamiltonian determines the variables, G and G, necessary 
for a complete quantum-mechanical description of this system. Hamilton's 
equations then yield 

rto = -2Ul + ±,-\^ 

tl A = - e 2 AG (15) 

or equivalently: 

A + e 2 GA = 

H§H(i) 2 -i^°' 

which correspond to (|TT|) and the expectation values of Eq. @ . 
The classical effective Lagrangian is 

This Lagrangian could also have been obtained using Dirac's action, 

r = J dt(*(t)\i^ - H\V(t)) = J dtL eS} (18) 

and a time-dependent Gaussian trial wave function as described in This 
variational method was used to study the quantum Henon-Heiles problem 
in a mean-field approximation 0. The Gaussian trial wave function is 
parametrized as follows 

tt(t) = [27rG(t)]- 1/4 exp[-(a: - ^t)) 2 ^" 1 ^)^ - iIl G (t)) + ip(t)(x - q(t))). 

We see that G(t) and nc(t) are the time dependent real and imaginary parts 
of the width of the wave function. One can prove for our problem that if the 
quantum oscillator starts at t = as a Gaussian, it is described at all times by 
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the above expression, where G(t) and Tlc(t) are totally determined by solving 
the effective Hamiltonian dynamics. (For our special initial conditions p(t) = 
q(t) = 0). Thus we find that our effective Hamiltonian totally determines the 
time evolution of the quantum oscillator. One interesting "classical" variable 
is the expectation value of the time dependent adiabatic number operator, 
which corresponds to the number of quanta in a situation where the classical 
A field is changing slowly (adiabatically). For the related field theory problem 
(see Section 3) of pair production of charged pairs by strong electric fields, 
this corresponds to the time dependent single particle distribution function 
of secondaries. To find the expression for the number of quanta, we begin 
with the wave function corresponding to a slowly varying classical field A: 

g(t) = exp [ - % J* u(t')dt'] / y/2u(t), 

in terms of which we can decompose the quantum operator via 

x(t)=g(t)b(t)+g*(t)tf(t). (19) 

Requiring the momentum operator to have the form 

p(t) = x(t) = git) b(t)+g*(t) b\t) 

by imposing g(t)b(t) + <7*(i)&t(i) = 0, and recognizing that g(t) and g*(t) 
satisfy the Wronskian condition by construction, then b(t) and tf(t) have 
the usual interpretation as creation and annihilation operators. That is, 
[x(t),p(t)] = i and [&(*), &t(i)] = 1. Also 

b(t)=i[g*(t)x(t)-g*(t)x(t)\. 

tf(t)b(t) can be interpreted as a time- dependent number operator for a slowly 
varying (adiabatic) classical field A. The time independent basis and the 
time dependent basis are both complete sets and are related by a unitary 
Bogoliubov transformation, b(t) = a(t) a + (3{t) a\ where 

a(t) = i[g*(t)f(t)-g*(t)f(t)\ 
Pit) = i[g*(t)f*(t)-g*(t)r(t)}, 

and where |ct(£)| 2 — |/9(t)| 2 = 1. If we choose for initial conditions, f2(0) = 
cj(0), Cl(0) = lj(0), then one finds that a(0) = 1 and (3(0) = 0. These are the 
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initial conditions appropriate to the field theory problem of pair production. 
The average value of the time-dependent occupation number is given by 



n(t) = (b\t)b(t)) = \f3(t)\ 2 = ^Qu)- 1 



• (20) 



Eq. ( p0|) allows us to compute the average occupation number of the 
system as a function of time. 





2 Solutions of the classical equations 

We first scale out the mass by letting t — > m~ l t, A —> wC 1 ! 2 A, G — ► vrT 1 G, 
and e — > em 3 / 2 . Then the scaled equations of motion are 

A + e 2 GA = 

2 i 

-^ + l + e 2 ^ 2 = 0. (21) 

In order to explore the degree of chaos as a function of (scaled) energy 
and coupling parameter e, we calculated surfaces of section and Lyapunov 
exponents. The surface of section is a slice through the three-dimensional 
energy shell |I| . That is, for a fixed energy and coupling parameter the points 
on the surface of section are generated as the trajectory pierces a fixed place 
(e.g. A = 0) in a fixed direction. The hallmark of regular motion is the 
cross section of a KAM torus which is seen as a closed curve in the surface 
of section. The hallmark of chaotic motion is the lack of any such pattern 
in the surface of section. In Fig. 1 we show a plot of a surface of section at 
E = 0.8 and e = 1. where regular and chaotic regions co-exist. 

The Lyapunov exponent provides a more quantitative, objective measure 
of the degree of chaos. The Lyapunov exponent, A, gives the rate of exponen- 
tial divergence of infinitesimally close trajectories ||. Although there are as 
many Lyapunov exponents as degrees of freedom, it is common to simply give 
the largest of these. For regular trajectories A = 0; for chaotic trajectories 
the exponent is positive. We define 

m = iim ^q+m-^) t (22) 

l*M \S\ 
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where z(zq, t) is a point in phase space at time t with initial position zq. Then 
the time evolution for ff(t) is 

?(t) = ff(t) ■ |^- 0j t) > (23) 

where 

l(t) = F(z(t),t) (24) 

are the full equations of motion for the system. The Lyapunov exponent is 
defined as 



A = lim - In 



ff(t) 



-7(0) ' (25) 

Appendix A of Ref. || provides an explicit algorithm for the calculation of 
all the Lyapunov exponents. Since we cannot carry out the t —>■ oo limit 
computationally, the regular trajectories are those for which A(i) decreases 
as while the chaotic trajectories give rise to A(t) that is roughly constant 
in time, as judged by a linear least-squares fit of log[A(t)] vs. log(t). 

We calculated the Lyapunov exponents for three values of the scaled cou- 
pling constant e (0.1, 1.0, 10.0) and for energies from 0.5 to 2.0. E = 0.5 
is the lowest energy possible, corresponding to the zero point energy of the 
oscillator; there is no upper limit on E. Fifty initial conditions were chosen 
at random for each energy bin of width 0.1 and coupling parameter. One 
relevant quantity to study is the chaotic volume, the fraction of initial con- 
ditions with positive definite Lyapunov exponents (corresponding to chaotic 
behavior). Errors in this quantity arise because of the finite number of ini- 
tial conditions chosen, and because the distinction between zero and positive 
exponents cannot be made with certainty at finite times. We found that for 
e = 0.1, more than 95% of trajectories were regular for all energies tested; 
for e = 1.0 and 10.0, there is a steadily increasing fraction of chaotic orbits 
between 0.5 < E < 1.25. For 1.25 < E < 2.0, more than 90% of these orbits 
are chaotic. 



3 Interpretation 

We may now ask what are the physical ramifications of our results. The 
system of equations studied here is the k = mode contribution to the 
problem of pair production of charged mesons by a strong electric field [^], 
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with E(t) = —A(t) being the value of the time evolving electric field. For 
that problem the equation for G(t) gets modified and becomes a function of 
the momentum k of the normal modes of the charged scalar field. Eq. (^) 
being replaced by oo 2 (t) — > u)\{t) = (k — eA(t)) 2 + m 2 . 
The semi-classical equation for A(t) becomes 

A = ej dk(k - eA(t))G k (t) , (26) 

which gets contributions from all modes. This system of equations is dis- 
cussed in detail in 0. When we sum over all the k modes, A(t) becomes 
a smooth function of time and is insensitive to initial data. However the 
number of particles produced in a narrow bin of momentum between k and 
k + dk depends only on Gk(t) and Gk{t). Only if one does a coarse grain- 
ing over momentum does one lose this sensitivity to initial data. Thus one 
should observe, as one counts the number of produced charged particles in a 
detector and increases the resolution, that the number of counts in a narrow 
momentum bin becomes a rapidly oscillating function of time whose behavior 
is chaotic. The expression for the number of particles in a given momentum 
bin k is 

n{k, t) = (An^r 1 ((Q k - co k ) 2 + 1 ^ - ^ j , (27) 

which should be compared with Eq. (0). The chaotic behavior of Eq. (p0|) 
is shown in Fig. 2. [Also see Fig. 3 of [[/]]. 
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FIGURE CAPTIONS 

Figure 1: A plot of the surface of section for energy = 0.8, e = 1.0, and 
A = 0. Each symbol represents a different trajectory. The one chaotic region 
is in the center of the plot. 

Figure 2: A plot of the occupation number given by Eq. (|20|) for energy = 1.8, 
e = 1.0, A(0) = 0, n G (0) = 0. The solid line is for for G(0) = 0.5; the dashed 
line is for G(0) = 0.5001. This plot shows the sensitivity to initial conditions. 
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